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Abstract 

We study the magnetic and transport properties of epitaxial graphene films in this letter. We 
predict enhanced signal of magnetic susceptibility and relate it to the intrinsic valley magnetic 
moments. There is also an anomalous contribution to the ordinary Hall effect, which is due to the 
valley dependent Berry phase or valley-orbit coupling. 
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Graphene, a monolayer carbon honeycomb lattice, has extraordinary electronic proper- 
ties[l,2]. Its conduction band and valence band form conically shaped valleys, touching at 
two nonequivalent corners of the hexagonal Brillouin zone called Dirac points. For optical 
and electronic application, it is desirable to make graphene semiconducting by opening a 
gap. A gap can be produced either by size quantization effects in graphene nano- ribbons [3], 
or by inversion symmetry breaking for epitaxial graphene grown on the top of crystals with 
matching lattices such as BN or SiC [4] . For the latter case, a significant energy gap has been 
predicted from ab initio calculation [5], although the experimental evidence for the gap is 
still under debate [6]. Experimentally, epitaxial graphene is relatively disordered and heavily 
doped because of close contact with the substrate. However, interesting physics can still 
manifest due to robust topological effects associated with the valley degree of freedom. 

In this Letter, we present our studies of magnetic and transport properties of epitaxial 
graphene systems. It is assumed that the interaction between the graphene layer and the 
substrate breaks the sublattice symmetry, which opens up a band gap at Dirac points [6]. 
We predict enhanced signals of magnetic susceptibility and anomalous contribution to the 
ordinary Hall effect. It is demonstrated that these extroadinary effects are due to the special 
valley degree of freedom in the system. Furthermore, we argue that both effects are robust 
against disorder, hence should be readily detected in experiment. 

Magnetic susceptibility. Recent studies on the orbital magnetism of 2D massless Dirac 
fermions have shown that the unique Landau-level structure gives rise to a singular diamag- 
netism[7]. The susceptibility becomes highly diamagnetic when the chemical potential is at 
band-touching point {/jl—O eV) where the density of states vanishes. 



Here, t is the nearest neighbor hopping energy, a is the lattice constant and (3 = l/Zc^T. T is 
the temperature and /cb is the Boltzman constant. It is still mysterious about the mechanism 
of such a singular behavior. A natural question is what happens to the susceptibility when 
a finite gap opens at the Dirac points. 

With a staggered sublattice potential, the low-energy effective Hamiltonian near the Dirac 
points is given by[8,9] 
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where Vq — y/3at/2 and H = k + eA is the reduced momentum operator in a magnetic 
field. A is the site energy difference between sublattices, which also corresponds to the band 
gap. Prom commutation relation [II-i-, nzp]==i=2ei?//?, where H-i- = U^^illy, we define ladder 
operators b and b'' such that: n_|_ = [y^ 2eB / h)b'^ and n_ = {^y2eB /fi)b. The wave function 
is constructed from the linear combination of eigenstates \n) of operator b. The Hamiltonian 
(2) is diagonalized to yield the Landau levels, which in conduction band is shown in Fig.l: 

, ,_/ r4 (n = 0), 

where = ±1 labels the two valleys. 

The magnetic susceptibility is defined as x = ~ (^|^) ' where the thermodynamic po- 
tential can be calculated from the formula F = — In [l + e^*^^""^"-*] . By using the 

n,Tz 

Euler-Maclaurin formula, we get the expand F as a power series with respect to B and the 
analytical expression of susceptibility is [10] 

^^^^ ~ AttH^A 1 + 2e/^'^Cosh(/3A/2) + e^/^M " ^ ■> 

Fig.2(a) shows the numerical result of the magnetic susceptibility at T = lOK. In the 
calculation we take A ~ 0.28 eV, t ~ 2.82 eV and a — 2.46A[5]. Compared with the zero 
gap case, it is observed that the large diamagnetic dip is still visible even when the gap 
opens, and it gets broadened in energy by the gap width. Disorder effects smooth out this 
curve but will not change the main features because susceptibility x is a thermodynamic 
property. We find that the integral of susceptibility over chemical potential, i.e. J 
is independent of both gap size and temperature. At zero temperature, the susceptibility 
becomes a square well shape, and vanishes when it is either electron- or hole-doped. 

The sudden jump of magnetic susceptibility at band edges signifies a large paramagnetism 
from the carriers. Indeed, if we calculate the magnetic susceptibility from Landau levels 
above the gap, then the conduction band contribution to x can be obtained as 

Prom the plot in Pig. 2(b), we can see that the large paramagnetic response from the conduc- 
tion electrons. What is the source of this paramagnetism? Note that our model does not take 
into account of spin degree of freedom, so it must come from the orbital motion of electrons. 



In fact, previous studies have shown that there is an intrinsic magnetic moment associated 
with the valley degree of freedom [11]. At conduction band bottom (ii — A/2), the moment 
equals the effective Bohr magneton /x^ = r^^^ with the effective mass m* = If one 

models this system by a 2D nonrelativistic electron gas with such a pseudo-spin degree of 
freedom, one would obtain both Pauli paramagnetic and Landau diamagnetic contributions 
recovering our expression Eq(5) above. 

Therefore, it is the valley magnetic moment that determines the magnitude of the sus- 
ceptibility signal. Using an experimentally measured energy gap of 0.28 eV for graphene on 
SiC[5], the effective mass is about 30 times smaller than the bare electron mass, hence the 
magnetic moment is about thirty times larger than the free electron spin magnetic moment. 

Semiclassical Landau levels. To further exemplify the role of the valley magnetic moment, 
we now closely examine the energy spectrum and analyze it via a semiclassical approach. 
If we place the exact Landau levels on the background of the original conduction band 
dispersion curve, one discovers that Landau levels from one valley start at the edge of the 
zero- field band (dashed curves in Fig.l), while those from the other valley start from one 
cyclotron energy above. The situation looks less odd if one shifts the bands by the Zeeman 
energy — m(k)-B, where m(k) = Tz 47-(^l^fp^272) is the valley magnetic momcnt[12]. Relative 
to the new bands (solid curves) , the Landau levels of both valleys now start at half of the 
cyclotron energy above the band edges, which is the famihar result for 2D free electron gas. 

Furthermore, one can obtain the Landau levels by semiclassically quantizing the cyclotron 
orbits in these modified bands. By taking into account the Berry phase correction, Onsager's 
quantization condition for the areas enclosed by the cyclotron orbits becomes vrfc^ = ^^y^[n-|- 
1 ~ ^2?][13]. In Rcf.[12], the Berry curvature of the conduction band at the two valleys has 
been calculated to be f2(k) — 2(A2+3fc°1^2t2)3/2 • "^^^ Berry phase r(A;) can then be obtained 
by integrating the Berry curvature over the area enclosed by the orbit. The energies of the 
modified bands at these quantized radius of the cyclotron orbits then yield the semiclassical 
Landau levels. Our result is shown in Fig. 3, which agrees very well with the exact Landau 
levels. However, we must emphasize that for the low lying Landau levels close to the band 
edge, the Berry phase effect is relatively less important than the Zeeman shift associated with 
the valley magnetic moments. The situation reverses at high energies, where the magnetic 
moment vanishes and the Berry phase approaches tt. 

Valley Polarization. The valley magnetic moment allows a direct coupling of the magnetic 
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field with the valley degree of freedom. We expect that an equilibrium population difference 
between the two valleys can be induced by applying a magnetic field. Analogous to the spin 
polarization concept, we define valley polarization to be Py = {n+ — n_)/ {n^ + n_), where 
n± is the densities of electrons in the valley with index = ±1. This can be calculated 
from the relative number of occupied Landau levels in the two valleys, resulting in the series 
of steps shown in Fig.4. Perfect valley polarization occurs when there is one Landau level 
occupied in one valley but none in the other(see Fig.l). 

We have also calculated the valley polarization semiclassically, in which the electron 
density from a given valley is obtained by integrating the Berry-curvature modified density 
of states up to the Fermi surface, i.e. (l + ^^^) [14]. The result is plotted as the 

black curves in Fig.4, and they smoothly go through the quantum steps. As expected, the 
valley polarization increases with the magnetic field (Fig.4 (a)), and the induced polarization 
is largest at the band edge where the valley magnetic moment is largest (Fig.4(b)). 

Anomalous Hall effect. It is now well established that in the presence of an in-plane 
electric field, an electron will acquire an anomalous velocity, proportional to the Berry 
curvature, in the transverse direction [13]. This leads to an intrinsic contribution to the 
Hall conductivity, a)^'* = 2(e^//l) J f {k)il{k) , where f{k) is the Fermi-Dirac distribu- 
tion function, and the factor of 2 comes from spin degeneracy. There is also a side-jump 
contribution [15] proportional to the Berry curvature at the Fermi surface when carriers 
scatter off impurities. Ignoring skew-scattering and other effects due to inter- valley scat- 
tering, a valley-dependent Hall conductivity at zero temperature is found in Ref.[12] as 
c_ff(Tz) = ^ ~ ~ • This result is independent of impurity density and 

strength. The valley dependence in the Hall current will lead to an accumulation of electrons 
on opposite sides of the sample with opposite valley indices. Clearly, if there is a net valley 
polarization, a charge Hall conductivity will appear upon the application of an electric field, 
(^xy — ~ 1^ <^A*- Here, p, is the chemical potential in the absence of the magnetic 

field, and S/i is the chemical potential difference between the two valleys. We can express 
S/j, in terms of the field-induced valley polarization as dfi = + n )^^B = jj^^^B, 

where n is the electron density and D{fi) is the density of states at the Fermi level. 

Assuming p^^ ^ p^y, the Hall resistivity may be expressed by —pxy = pix'^xy = JahB, 
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where ^ah is the anomalous Hall coefficient [16] given by 

with me being the bare electron mass. Figure 5 shows the comparison of the anomalous Hall 
coefficient with the ordinary one ^qh {—~^) zero temperature. It is shown that with 
increase of chemical potential, ^ah changes from positive to negative and has a curvature 
opposite to that of ^oh- Moreover, we find that ^ah is comparable with ^qh for a typical 
value of p-r-r=10 A;J1[17] and becomes dominant for larger resistivity. It indicates that Hall 
coefficient may be detected as a signal of Berry curvature experimentally. 

Valley- orbit coupling. The traditional theories of anomalous Hall effects are all based on 
two ingredients: spin polarization and spin-orbit coupling [18]. In our system, spin-orbit 
coupling is extremely weak [19], so we have ignored it from the very beginning. Because 
of the magnetic moment associated with the valleys, we are able to produce a population 
imbalance between the two valleys, i.e., a valley polarization, by a magnetic field. Is there 
a valley-orbit coupling that underhes the anomalous Hall effect discussed in this work? 

From a phenomenological point of view, there is indeed a valley-orbit coupling, because 
electrons in two valleys do have opposite anomalous velocities under an electric field. In fact, 
the analogy with spin-orbit coupling can be made more precise if we consider the effective 
Hamiltonian for the conduction band. Ref.[13] has provided a recipe to construct the effective 
Hamiltonian from three basic ingredients: the band energy, the magnetic moment, and the 
Berry curvature. The magnetic moment gives rise to the Zeeman coupling between the 
magnetic field and the valley index. There is also a dipole-like term proportional to the 
electric field eE • R, where R is the Berry connection potential defined by x R = 
Using the expression of the Berry curvature for the graphene system given by Eq.(2), we 
find that R = — ^ ^1 — -^^5=5=5^^ k x z. Therefore, the low-energy effective Hamiltonian 
[20] can be written as 

T en? 

H^n = £o(k) - r,/.^B • z - • (k X E), (6) 

where c* = ^^^^ is the electron speed at Dirac point of the gapless graphene. The third 
term is the valley-orbit coupling, which has the similar form as the spin-orbit coupling. 
The effective Hamiltonian resembles closely the Pauli Hamiltonian for free electrons in the 
non-relativistic limit. 
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In summary, we have studied the magnetic susceptibihty and anomalous Hall effect in 
epitaxial graphene. Our results show that a large diamagnetism will result from the large 
valley magnetic moment of graphene electrons. We also predict an anomalous Hall coeffi- 
cient, which can dominant that of the ordinary Hall effect, due to a valley dependent Berry 
phase or valley-orbit coupling. 
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FIG.l: (Color online). The Landau levels from exact quantum calculations for Tz — +1 

(a) and Tz — —1 (b) valleys. For comparison, the energy dispersions £o(k) (dashed hne) and 
£:(k) — £o(k) — m(k) ■ B (solid curves) are shown. Parameters are A= 0.28 eV , t— 2.82 eV 
and lattice constant a=2.46 A. 

FIG. 2: (Color online), (a) Magnetic susceptibility as a function of chemical potential. 

(b) Magnetic susceptibility from conduction electrons only. Here, T — 10 K and A= 0.28 
eV. 

FIG. 3: (Color onhne). The field dependence of semiclassical Landau levels for Tz — +1 
(red circle) and — —1 (blue triangle) valleys. For comparison, the exact quantum Landau 
levels are also shown (solid line). 

FIG. 4: (Color online). (a)The variation of valley polarization Py with the magnetic 
field B {/I— 0.2 eV). (b)The variation of valley polarization Py with the chemical potential /i 
(B=0.2 T). Both the quantum calculations (red) and semiclassical results (black) are shown. 

FIG. 5: (Color onhne). The variation of anomalous Hall coefficient ^ah, ordinary Hall 
coefficient ^yoH and the total {'Jah+^oh) with the chemical potential /i. Here, p^x — 10 kD,, 
A = 0.28 eV, t = 2.82 eV and a = 2.46A. 
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FIG. 1: (Color online). The Landau levels from exact quantum calculations for = +1 (a) 
and Tz = —1 (b) valleys. For comparison, the energy dispersions eo(k) (dashed line) and e(k) = 
eo(k) — m(k) • B (solid curves) are shown. Parameters are A= 0.28 eV, t= 2.82 eV and lattice 
constant a=2.46 A. 
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FIG. 2: (Color online), (a) Magnetic susceptibility as a function of chemical potential, (b) Magnetic 
susceptibility from conduction electrons only. Here, T = K and A= 0.28 eV . 
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FIG. 3: (Color online). The field dependence of semiclassical Landau levels for = +1 (red circle) 
and Tz = —1 (blue triangle) valleys. For comparison, the exact quantum Landau levels are also 
shown (solid line). 
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FIG. 4: (Color online). (a)The variation of valley polarization Py with the magnetic field B {fj,=0.2 
eV). (b)The variation of valley polarization Py with the chemical potential fj, (B=0.2 T). Both the 
quantum calculations (red) and semiclassical results (black) are shown. 
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FIG. 5: (Color online). The variation of anomalous Hall coefficient "fAH-, ordinary Hall coefficient 
^OH and the total {'^ah+Ioh) with the chemical potential //. Here, pxx = 10 fcfi, A = 0.28 eV , 
t = 2.82 eV and a = 2.46A. 
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